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Abstract 

The Casimir force on a D-dimensional sphere due to the confinement of a 

massless scalar field is computed as a function of D, where D is a continuous 

variable that ranges from — oo to oo. The dependence of the force on the 

dimension is obtained using a simple and straightforward Green's function 

technique. We find that the Casimir force vanishes as D — > +oo (D non-even 

integer) and also vanishes when D is a negative even integer. The force has 

simple poles at positive even integer values of D. 
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I. INTRODUCTION 



In recent papers it was proposed that the dimension of space-time could be used 
as a perturbation parameter in quantum field theory calculations. The advantage of such 
an approach is that analytical (nonnumerical) results can be obtained which are nonpertur- 
bative in the coupling constant. This procedure was used to obtain the Green's functions 
of self- interacting scalar quantum field theory in the Ising limit |5]||. One can also perform 
dimensional expansions in inverse powers of the dimension. Such expansions have proved 
useful in atomic physics calculations 0. 

These perturbative investigations have led to and motivated analyses of the dimensional 
dependence of various physical systems. Such investigations are useful because by identifying 
the singularities in the complex-dimension plane one can predict the radius of convergence 
of a dimensional expansion. The dimensional dependence of some elementary quantum- 
mechanical and field-theoretic models is described in The dimensional dependence of 
classical physical models has also been investigated; for example, in |||| the dimensional 
dependence of probabilities in models of random walks was elucidated. 

In this paper we investigate the dimensional dependence of the Casimir force per unit 
area, F/A, on a spherical shell of radius a in D space dimensions. Specifically, we study the 
Casimir force that is due to quantum fluctuations of a free massless scalar field satisfying 
Dirichlet boundary conditions on the shell. 

An interesting investigation of the dependence of the Casimir force per unit area upon 
the spatial dimension is already in the literature [10|. Ambj0rn and Wolfram examined 
the case of infinite parallel plates embedded in a D- dimensional space and separated by a 
distance 2a; that is, there is one longitudinal dimension and D — 1 transverse dimensions. 
Their result is 

F/A = -a- D - x 2- 2D -\-^l 2 D T (^Y^j CP + 1), (1.1) 

which we have plotted in Fig. 1. Note that F/A has a simple pole (due to the Gamma 
function) at D = — 1. However, F/A is not infinite at the other poles of the Gamma function, 



which are located at all the negative odd integral values of D because the Riemann zeta 
function vanishes at all negative even values of its argument. One interesting and well-known 
special case of ( |1 . 1| ) is D = 1: 

F/4~ - ~ (1-2) 

where the negative sign indicates that the force is attractive. We mention this case here 
because the spherical geometry considered in the present paper coincides with the slab 
geometry of [|l(]] when D = 1; we recover the result (|1.2| ) as a special case in Sec. [II . 



This paper is organized very simply. In Sec. |IJ we review the Green's function formalism 
required to obtain the Casimir force. Then we apply this formalism in D- dimensional space 
to obtain an expression for the Casimir force per unit area on a D-dimensional spherical shell. 
This expression takes the form of an infinite sum of integrals over modified Bessel functions; 
the dimension D appears in the orders of the Bessel functions. In Sec. [TTT| we examine this 
expression for the Casimir force per unit area in detail. We show that each term in the series 
exists (each of the integrals converges) and we show how to evaluate the sum of the series 
numerically for all real D. When D > the Casimir force per unit area is real; the force is 
finite except when D is an even integer. When D < the Casimir force is complex; there 
are logarithmic singularities in the complex-D plane at D — 0, — 2, — 4, — 6, .... 

II. MATHEMATICAL FORMALISM 

The calculation in this paper of the Casimir force on a spherical shell relies on the use of 
Green's functions to represent vacuum expectation values of time-ordered products of fields. 
The Green's function is used to obtain the vacuum expectation value of the stress-energy 
tensor, from which we will derive the Casimir force. The formalism used here was developed 



in 



11|-[T3[]. We summarize the formalism below. 

A free massless scalar field y?(x, t) satisfies the Klein-Gordon equation 

d 2 



()f2 V^(x,t) = 0, (2.1) 



where x is a D-dimensional position vector. (Initially, we will think of D as a positive integer; 
however, once we have derived the radial equation for the Green's function we will be able 
to regard the parameter D as a continuous variable.) The quantum nature of the Casimir 
force arises from the constraint that <p(x,t) satisfies equal-time commutation relations: 

b(x,t),0(x',t)]=^)(x-x'). (2.2) 

The two-point Green's function G(x, t; x', t') is defined as the vacuum expectation value 
of the time-ordered product of two fields: 

G(x,t;x',t') = -i(Q\T<p(x,t)<p(y,t')\Q). (2.3) 

By virtue of (|2~T| ) and ( |2.2| ), the Green's function G(x, t; x', t') satisfies the inhomogeneous 
Klein-Gordon equation 

(j^ - V 2 ) G(x, t; x', t') = -<^(x - x')5(i - t'). (2.4) 

We will solve the above Green's function equation by dividing space into two regions, region 
I, the interior of a sphere of radius a and region II, the exterior of the sphere. On the sphere 
we will impose Dirichlet boundary conditions 

G(x,t;x',t') | |x|=a =0. (2.5) 

In addition, in region I we will require that G be finite at the origin x = and in region II 
we will require that G satisfy out going- wave boundary conditions at |x| = oo. 
The stress-energy tensor T Ml/ (x, t) is defined as [14[ 

T^(x, t) ee d*V(x, t)<9V(x, t) - l -g^d x y{^ t)dV(x, t). (2.6) 

The radial Casimir force per unit area F/A on the sphere is obtained from the radial-radial 



component of the vacuum expectation value of the stress-energy tensor [11]: 



FM=(o|r--r-|o)| r=a . (2.7) 
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To calculate F/A we exploit the connection between the vacuum expectation value of the 
stress-energy tensor T^(x, t) and the Green's function at equal times G(x, t; x', t): 



F/A 





-G(x, t; x, t) m - — — G(x, t; x', t\ 



Or Or' 



(2.8) 



x=x', |x|=a 

To evaluate the expression in ( |2.8| ) it is necessary to solve the Green's function equation 
^|). We begin by taking the time Fourier transform of G: 

(2.9) 



/oo 
d t e -Mt-t') G ^t;x.',l/). 
-oo 

The differential equation satisfied by C? w is 

(^ + V 2 )^(x;x') = 5 (D) (x-x'). 



(2.10) 



To solve this equation we introduce polar coordinates and seek a solution that has cylin- 
drical symmetry; i.e., we seek a solution that is a function only of the two variables r = |x| 
and 9, the angle between x and x' so that x • x' = rr' cos 9. In terms of these polar variables 
( P-lOj ) becomes 



2 D - 1 sin 2 ~ D 9 



*~ Or 2 ^ r Or 







<f(r - r')5(^r (^i) 

27r (D-l)/2 r D-l sin D-2^- 

(2.11) 



Note that the D-dimensional delta function on the right side of ( |2.10|) has been replaced by 
a cylindrically-symmetric delta function having the property that its volume integral in D 
dimensional space is unity. The D-dimensional volume integral of a cylindrically-symmetric 
function f(r, 9) is 

2tt{D-1)/2 poo pit 

r(D n / drr - 1 d9 sin D " 2 9 f(r, 9). (2.12) 
1 {.—) Jo Jo 



We solve ( |2.11 ) using the method of separation of variables. Let 



g u ( r y,e) = A(r)B( z ) 



(2.13) 



where z = cos9. The equation satisfied by B(z) is then 



'1 _ Z *)S^ - Z [B _ i) A. + n („ + d - 2) 



BW=0, 



(2.14) 



dz 2 ' ' <fc 

where we have anticipated a convenient form for the separation constant. The equation 
satisfied by A(r) is 

d 2 D - 1 d n(n + D - 2) 



A(r)=0 (r^r'). 



(2.15) 



<ir 2 r dr r 2 
The solution to fl2.14|) that is regular at \z\ = 1 is the ultraspherical (Gegenbauer) polynomial 



in 



B{z)=C^ 1+D ' 2 \z) (n = 0, 1,2,3, ...). 



(2.16) 



The solution in region I to (2.15) that is regular at r = involves the Bessel function 16 



A(r) = r 1 - fl / a J B _ 1+ a(|o;|r). 



(2-17) 



In ( p,17| ) we assume that D > 2 in order to eliminate the linearly independent solution 
A{r) = r 1 ~ D / 2 Y n _ 1+ D (\u\r), which is singular at r = for all n. The solution in region II 
to ( |2.15| ) that corresponds to an outgoing wave at r = oo involves a Hankel function of the 
first kind |§ 



A(r)=r 1 - D / 2 H^ 1+R (\uj\r). 



(2.18) 



The general solution to ( 2.11]) is an arbitrary linear combination of separated-variable 
solutions; in region I the Green's function has the form 



Q u {rS,0) = Y,^r 1 - D/2 J n _ 1+§ (\uj\r)Ci- 1+D / 2 \z) (r < r' < a) 



(2.19a) 



n=0 



and 



n=Q 



g u (r,r>,6) = Y j r l ~ D > 2 [b n J n _ 1+ »(\u;\r) + c n J_ n+1 _ f (Mr)] Ct 1+D/2 \z) (r> < r < a). 

(2.19b) 



[Note that J v {x) and J- V (x) are linearly independent so long as v is not an integer. Thus, 
( |2.19b| ) assumes explicitly that D is not an even integer.] The general solution to fl2.11| ) in 
region II has the form 
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g u (r,r',e) = f^dy- D / 2 H^ +R (\u\r)Ct 1+D/2) (z) (r > r' > a) (2.20a) 



n=0 



and 



^(r,^e) = ^r 1 ^/ 2 le^ 1+ ^(|^|r) + / n ^ i+ ^(|^|r) C^+^C*) 



(2) 



n=0 



(r' > r > a). 
(2.20b) 



The arbitrary coefficients a n , b n , c n , d n , e n , and f n are uniquely determined by six 
conditions; namely, the Dirichlet boundary condition (|2.5| ) at r = a, 

(2.21a) 



b n J n - l+ D (\u\a) + c n J_ n+l _D{\u)\a) = 



and 



e n H^ (\u\a) + f n H^>{\u\a) = 



(2) 



the condition of continuity at r = r', 



a n J n _ 1+ D_(\uj\r') = b n J n _ 1+ D(\u\r') + c n J„ n+1 n(\u\r') 



and 



r(i) 



(2) 



(2.21b) 



(2.21c) 



and the jump condition in the first derivative of the Green's function at r = r', 



(2.21d) 



b n J' n _ 1+ „(\u\r') +c n J'_ n+1 _ R (\u\r') - a n J' n _ 1+!1 {\u\r>) = {2n + D ./f^ 2 ^ (2.21,0 



4(7rr') 2 | co' | 



and 



e n ^' 1+c (|.|r') + / n iff i+D ( M r') - d^' (Mr') = ^ + ^ ■ (2.21f) 

n i+ 2 n i+ 2 n i+ 2 4(7rr') 2" | a; | 

Here we have used the orthogonality property of the ultraspherical polynomials |15 

• l 
-i 



dz{l-z 2 Y-y 2 C^{z)C%\z) = 2l ~u* V{ \t^ 5n m (a^O), (2.22) 



n\ {n + a)T 2 (a) 



the value of the ultraspherical polynomials at z = 1, 

^W = ^W (a#0) ' (2 ' 23) 

and the duplication formula, F(2x) = 2 2x ~ 1 T(x)T(x + 1/2) /y^. 

Having determined the coefficients in the expressions for the Green's function, we can 
immediately evaluate the right side of (|2.8|) . The contribution to F/A from the interior 
region (region I) is 



{F/A), 



™ {n-l+J )Y(n + D-2) 



n=0 



2^7r^a^n!r(^) 



da; 



Ma 



hi 

J n _ 1+ D{\u\a) 2 



. (2.24a) 



The contribution to F/A from the exterior region (region II) is 



(WIA\ ,^ (n-l + f)r(n + i9-2) 



do; 



a; 



ai/^' (Ma) 



H 



(i) 



, (|a>|a) 



D 



. (2.24b) 



n=0 * « - u, ii! r ( i ^) 

The integrals in (|2.24 ) are oscillatory and therefore very difficult to evaluate numerically. 
Thus, it is advantageous to perform a rotation of 90 degrees in the complex-o; plane. The 
resulting final expression for F/A is 



p/A ^ (n-l + §)r(n + D-2) 



n=0 



% 2 a D+1 n 



dx 



xI ' n _ l+ o(x) 
l n . 1+ o(x) 



xK' 



+ 



X 



K 



71-1 + 



D (x) 



+ 2-D 



(2.25) 



The D = 2 result, where the n = term appears with weight ~, was derived in This 
result can be recovered by setting D = 2 + e and letting e tend to 0. 



III. NUMERICAL EVALUATION OF F/A 



Our objective now is to evaluate the formal expression in ( |2.25| ) for arbitrary dimension 
D. Recall that this expression was derived under the assumption that D > 2 and that D is 



not an even integer. However, we will now seek an interpretation of ( |2.25| ) that is generally 
valid; to do so we will apply a summation procedure that enables us to continue fl2.25| ) to 
all values of D. 
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The expression in ( 2.25Q does not exist a priori for all D. Furthermore, as we will see, 
the individual terms in the series, which are integrals in x, do not exist. Fortunately, it is 
possible to modify the terms in the series so that the integrals do exist; this modification 
requires a delicate and detailed argument. However, there is one simple case, namely that 
for which D = 1, where the series fl2.25p is well-defined and easy to evaluate. We examine 
this case in the next subsection. 



A. Special case D = 1 



When D = 1 the series in Q2.25| ) truncates after two terms. This happens because of the 
identity 



r(n + L>-2) 1 1 
When this identity is inserted into the sum in ( |2.25| ) we obtain 



F/A\ 



1 



D=l 



Ana 2 j 
Next we use the identities 



dx 



xl[ (x) xl' i (x) xK[ (x) xK' ! (x) 



+ 



+ 



h(x) I_i(x) Ki(x) 

2 2 2 



+ 



K_i(x) 



+ 2 



(3.1) 



(3.2) 



sinhx coshx 



, Ki{x) = Ki(x) = 
x 2 2 Jx 



(3.3) 



to reduce (|3.2|) to the (convergent) integral 

1 



F/A\ 



D=l 



Ana 2 



dx 



d j / coshxsinhx . } 

dx \ x 2 e 2x 



Ana 2 
1 



If 00 d 
m 2 dXX dx' ln ( 1 - e ~^ 



dx- 



x 



7TO j 
7T 



96a 2 



(3.4) 



This result agrees with the well-known result given in (|1.2j ). 
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B. Convergent reformulation of ( |2.25| ) 

In this subsection we modify the form of the series ( |2.25| ) so that each term in the series 
exists and we apply a summation procedure to evaluate the the resulting series numerically. 
We begin by rewriting ( f2.25[ ) in the form 

F/A = - Y {U ~ 1 tP T{n + D ~ 2) rdx\x±\n(l n l+ »(x)K n l+ o(x))+2-D 

(3.5) 

In this form it is easy to investigate the convergence of the individual integrals in the series. 
To do so we recall the asymptotic behavior as x — > oo: 

I u {x)K v (x) ~ — (x -» +00). (3.6) 



From Q3.6|) it is clear that the integrals in ( |3.5|) do not converge except for the special case 
D — 1. However, as we will now argue, one can replace the quantity 2 — D in ( p.5|) by 1 
without changing the value of F/A, provided that D < 1. This replacement will render the 
integrals convergent. 
Consider the series 



oc 



V r(n + a) - (3-7) 

n=0 

This series converges so long as a < and a 7^ —N, where N — 1, 2, 3, ... (so that 
individual terms in the series exist). This series can be summed in closed form because it is 
a special case of the binomial expansion 

n=0 U ' 

Note that if we let x — > 1 we obtain the identity 

T(n + a) , 

n=0 



One can also show that the identity ( ^.9| ) holds in the limit as a; approaches a negative 
integer — N. To do so we let a = —N + e, where iV = 1, 2, 3, .... We then decompose the 
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series into two parts, the first whose terms are finite and the second whose terms diverge as 
e -> 0: 

~ n \ i N (-l) n 

^ (n + N+ 1)! + 7^ (N-n)\T(l + n-ey 

n=0 n=0 ' 

The first sum in (|3.10|) can be easily evaluated as 

1 



NNl 



(3.11) 



If the second sum in fl3.10|) is expanded in a Laurent series in e, the coefficient of 1/e vanishes. 



However, the coefficient of e° is given by 



n=0 

Finally, using the integral representation 



1 » / A/A 

MB" 1 )" J *(»+!)■ < 3 ' 12 > 



-1 i j.n 

V»(n + l) = -7+ / -, (3.13) 



we show that (|3.12|) cancels ( p. 11 



This argument shows that if the integral in ( |3.5| ) is replaced by 1 then the sum vanishes: 

2 -^-4 ^ = (D<1). 3.14 

/ — ' n! 

n=0 

It follows that in the region D < 1 (if we sum first over n) we can add any constant to 
the integrand in each term in the series ( |3.5| ) without changing the value of the sum. We 
conclude that we may replace 2 — D by 1 in ( |3 . 5| ) . Our new improved expression for the 
Casimir force per unit area is thus 

f ' a = t ^rfi^rif [ dx ^K-fW^-fW) ■ < 3 ' i5) 

For D > each term in this series exists for all n. (For D < there is yet another subtlety 
that we will address shortly.) Before we proceed, we must emphasize that while ( |3.15| ) has a 



different and more compact form than that in ( 2.25|) we have not changed the value of F/A; 



we have in effect added zero to the series representing F/A. 
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Unfortunately, the formula in (|3.15|) is still not satisfactory because the series does not 
converge. To examine the convergence of this series we need to know the asymptotic behavior 
of the integrals for large n. We make use of the uniform asymptotic approximation to the 
product I v {vx)K v (i>x): 

I v (ux)K„(ux) ~ — (l + — 2 + ...J (i/->oo), (3.16) 

where £ = (1 + x 2 ) -1 / 2 . This asymptotic behavior implies that the integral, which we will 
abbreviate by Q n , in the nth term in ( |3.15| ) grows linearly with increasing n: 

/oo f v \ 35 \ 

dx \n[2xI v {x)K u {x)\ ~ 7T I - + — - 3276gu3 +■■■) {n -> oo), (3.17) 

where z/ = n — 1 + y . Because of this linear growth in n it is apparent that the series in (p. 15 ) 
does not converge if D > except for the special case D = 1, where the series truncates! 

To solve this problem we introduce an analytic summation procedure based on the prop- 
erties of the Riemann zeta function. Specifically, we consider the leading large-n behavior 
of the summand in Q3.15 ): 

(n - 1 + f )r(n + D - 2) 1 D _ x 



' 2 -Q n 5-5 (n u - 1 + ...) (n->oo). (3.18) 

Then, regarding the parameter .D as being less than 0, we sum the expression on the right 
side of ( |3.18j ) over n from 1 to oo. This gives 

5-i — C(l-D), (3.19) 

which is a well-defined function of D. We now add Q3.19 ) to ( 3.15|) and correspondingly 
subtract the right side of fl3.18| ) from each term (except the n = term) in the expression 
for F/A in ( |3.15| ). This produces a new series for F/A that is convergent for D < 1: 

D 



F/A = TrrS^r/ dx\n[2xl ,,d(x)K , + D (x 



n=l 



7T D . n-l + f T(n + L>-2 

-n^' 1 + 

2 n! 



x J™ dx\n(2xI n _ 1+ n(x)K n _ 1+ n(x)) - 2gr /^ix C(l--P)}- (3.20) 
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We have finally achieved our objective; we have obtained a convergent series representation 
for F/A for a finite range of positive D, namely, < D < 1. Each term in this series exists. 
In the next subsection we will use this series to calculate analytically the Casimir force on 
a sphere in zero dimensions. 



C. Casimir force on a zero-dimensional sphere 



If we substitute D = in ( |3.20|) we find that F/A = oo; this is because the Riemann 



zeta function ((1 — D) is singular when D — 0. The divergence in F/A is a consequence of 
that fact that the surface area of a zero- dimensional sphere is zero. Hence, we will compute 
the Casimir force, rather than the Casimir force per unit area. The equation for F is 
obtained by multiplying (|3.20|) by the surface area of a D-dimensional sphere of radius a, 
A = 2a D - 1 7r D / 2 /r( J D/2), 



F =^{1 dxln(2 X I_ 1+f (x)K_ 1+ „ 
2 ~ 



x 



k D _ x , [n-l + §)V(n + D-2) 



-n + 

2 //! 

71 



dx In (2xI n ^ 1+ D(x)K n _ 1+ D{x 



r(D-i) 



((1-D)\. (3.21) 



We can now let D tend to in ( |3.21|) . We obtain the result 

1 



\D=0 



(3.22) 



2a 2 ' 

where we have used ((z) ~ -^r as z — > 1. 

Note that we could not have obtained this result from ( |3.15| ). Indeed, if we naively let 
D — > 0+ in the formula obtained by multiplying (|3.15| ) by the surface area A of the sphere 
we appear to get the value 0. This is because only the n = and n = 2 terms survive 
in this limit and these two terms cancel as a result of the identities I- n {z) = I n {z) and 
K- U (z) = K v (z). However, the result F = at D = is incorrect because the series in 
( P-15| ) does not converge. 
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D. Numerical results 



The expression in (|3.20| ) may in principle be used to compute F/A numerically for D < 1; 
to wit, we may evaluate the integrals for a large number N of terms in the series, compute 
the iVth partial sum, and extrapolate the result to its value at iV = oo. However, this 
procedure is rather inefficient because the sum in ( |3.20|) is very slowly converging. Therefore, 
to prepare for evaluating F/A we subtract not just the one term in ( |3.18| ) but many terms 
in this asymptotic expansion. The first three terms are 



n-l + §)T(n + D-2) 1 r D _ x p- !)(£>- 2) 



2AD 4 - 176L> 3 + 504D 2 - 688L> + 387 n , 

-n u ~ 6 + . 



n H n 



D-2 



- m (n-oo). (3.23) 

If we use K terms in this asymptotic expansion we then have K corresponding Riemann 
zeta functions appearing in the final form of the series. The series converges more rapidly 
(the nth term in the series vanishes like n D ~ K ~ l ) and it also converges for a larger range of 
the dimension: D < K. We have used this method to graph F/A and F as functions of D 
(see Figs. 2 and 3). 

From Figs. 2 and 3 it appears that F and F/A are singular at D = 2 and D — 4. If fact, 
as we will now explain, the Casimir force is singular at all even positive integer values of D; 
F and F/A have simple poles at D = 2N, N = 1, 2, 3, .... To verify this, we examine the 
generalization of ( |3.20|) obtained by making many subtractions of the asymptotic behavior 
in (|3.23|) . This formula for F/A will contain many Riemann zeta functions, one for each 
subtraction. The kth zeta function will have the form ((k — D). Furthermore, if k is even 
the coefficient of ((k — D) contains the factor (k — 1 — D). (No such factor occurs if k is 
odd.) Thus, when D = k — 1 and k is an even positive integer the simple pole of the zeta 
function is cancelled by this factor and when D = k — 1 and k is an odd positive integer, 
the simple pole persists. 

As explained above, the Casimir force is finite at all odd-integer dimensions. For example, 

^' D=3 = i-0.0028168..., (3.24) 



a- 
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where the positive value indicates that the force is repulsive (tends to inflate the sphere) . The 
numerical value in ( |3.24| ) is much smaller than that obtained by Boyer [[H]] for the case of an 
electromagnetic field confined in a three-dimensional spherical cavity (F = 0.046176 . . . a~ 2 ) 
and that obtained by Milton for the case of a spinor field confined in a three-dimensional 
spherical cavity (F = 0.0204 . . . a~ 2 ). 

Although we have not proved it, it does appear from Fig. 3 that F — > as D — > oo. This 
is probably associated with the fact that the volume and surface area of a D-dimensional 
sphere of radius a tend to as D tends to oo. 



E. Casimir force for negative dimension 



For all odd-integer D < 1 the series in ( |3.15 ) truncates and thus it is not necessary to 



subtract off the large-n behavior. This truncation occurs because of the identity 

which is the generalization of ( |3.1| ) . If we use this identity at iV = 1 we obtain the following 
integral representation for F/A at D — — 1: 



1 r°° 

F/A\ D= ^ = -- I dx In 



o 



1 \ 2 / sinh x \ ( cosh x 
1 + - 1- 



(3.26) 



x ) \ x cosh x J \ x sinh x / 

where we have inserted the expressions for the half-odd integer modified Bessel functions. 
An interesting aspect of this integral representation for F/A is that the argument of the 
logarithm has a zero for a positive real value of x. This zero may be traced to the positive 
zero of the function /_ 3 / 2 (a;). Thus, F/A\b=-i is complex! Because the contour of integration 
passes under the zero, we find here F/A\d=-i = 0.65382 + il. 88445. 

In general, for all D < 0, D ^ —2N with N = 1, 2, 3, . . ., the argument of the logarithm 
in the integrand of ( p. 15] ) always has a zero. Hence, the analytic continuation of F/A to 
negative values of D is complex. The zero of the argument of the logarithm comes about 
because I u (x) has a real positive zero when —2m < v < —2m + 1, where m is a positive 
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integer. Only a finite number of integrals in the series in ( |3.15| ) are complex. However, 
as D becomes more negative, there are more and more complex integrals in the series. In 
particular, each time D decreases past a negative even integer one additional integral in 
the series ( |3.15| ) becomes complex. Thus, in the complex- .D plane, F/A has branch cuts 
emanating from the points D = —2k, k = 0, 1, 2, 3, . . .. 

It is remarkable, however, that exactly at the negative even integers, it is possible to 
evaluate the Casimir force F; we find that at these points F = 0. This is because the series 
( 3.15Q truncates after a finite number of terms for these values of D, and the remaining 
terms cancel in pairs. In Fig. 4 we plot KeF for —5 < D < 5. Figure 4 illustrates one 
interesting aspect of the Casimir force, namely, the erratic fluctuations in the sign of F. 
The sign of the Casimir force is extremely difficult to understand intuitively — we know of 
no simple physical argument that predicts whether the force is attractive or repulsive. 

We thank the US Department of Energy for financial support. 
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FIG. 2. A plot of the Casimir force per unit area F/A for < D < 5 on a spherical shell. 
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FIG. 3. A plot of the Casimir force F for < D < 5 on a spherical shell. 
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